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ABSTRACT 

An  elastic-viscoplastic  model  for  describing  the  anisotropic  high-strain  rate  behavior  of 
both  low-symmetry  and  high-symmetry  textured  materials  is  proposed.  Yielding  is 
described  using  a  recently  developed  criterion  which  can  capture  simultaneously 
anisotropy  and  compression-tension  asymmetry  associated  with  deformation  twinning. 
The  anisotropy  coefficients  as  well  as  the  size  of  the  elastic  domain  are  considered  to  be 
functions  of  the  accumulated  plastic  strain.  The  specific  expressions  for  the  evolution 
laws  are  determined  using  a  multi-scale  methodology,  i.e.  experimental  measurements  of 
crystallographic  texture  and  uniaxial  stress-strain  curves,  polycrystalline  calculations,  and 
macroscopic  scale  interpolation  techniques.  An  overstress  approach  is  used  to  incorporate 
rate  effects  in  the  formulation.  Applications  of  the  model  to  the  description  of  the  high 
strain-rate  response  of  low-symmetry  (hexagonal-close-packed  zirconium)  and  high- 
symmetry  (body-centered-cubic  tantalum)  pre-textured  metals  are  presented.  The  very 
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good  agreement  between  the  simulated  and  experimental  post-test  geometries  of  the 
Taylor  impact  specimens  in  terms  of  major  and  minor  side  profiles  and  impact-interface 
footprints  shows  the  ability  of  the  model  to  describe  with  fidelity  the  differences  in  the 
evolution  of  anisotropy  between  zirconium  and  tantalum. 


Keywords:  anisotropic  material  (B),  elastic-viscoplastic  matetrial  (B),  crystal  plasticity 
(B),  impact  testing  (C),  finite  elements  (C) 
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1.  INTRODUCTION 


The  importance  of  an  adequate  description  of  plastic  anisotropy  has  been  demonstrated  in 
many  low  strain-rate  forming  applications  (e.g.  Cazacu  and  Barlat,  2003;  Yoon  et  al., 
2006;  Soare  et  al.,  2006,  etc.).  Recent  efforts  in  the  development  of  computational  models 
to  describe  temperature  and  strain-rate  effects  on  the  inelastic  response  have  given  rise  to 
robust  predictive  methods  for  simulating  the  anisotropic  high  strain-rate  behavior  of 
body-centered-cubic  (bcc)  and  face-centered-cubic  (fee)  polycrystals.  Maudlin  et  al. 
(1999a,b)  employed  an  anisotropic  representation  of  the  yield  surface  obtained  from 
polycrystal  calculations  to  model  the  mechanical  response  of  a  pre-textured  bcc  tantalum 
sheet.  Anand  et  al.  (1997)  developed  finite-element  (FE)  polycrystal  models  for 
describing  the  high  strain-rate  and  large  deformation  behavior  of  both  bcc  and  fee 
polycrystals.  These  models  describe  qualitatively  well  the  important  mushrooming  of  the 
Taylor  impact  specimens.  In  contrast,  the  prediction  of  the  geometry  changes  observed  in 
hexagonal-close-packed  (hep)  specimens  remains  a  challenge  (Maudlin  et  al,  1999b). 
Metallographic  investigations  have  shown  that  while  deformation  in  hep  materials  under 
quasi-static  and  dynamic  loading  includes  both  slip  and  deformation  twinning,  the 
propensity  of  twinning  is  higher  for  high-rate  loading  and/or  low-temperature 
deformation.  Therefore,  in  order  to  describe  with  accuracy  the  dynamic  response  of  hep 
materials  it  is  imperative  to  account  for  the  most  important  sources  of  anisotropy  in  the 
given  material:  slip  and/or  twinning  activity,  substructure  evolution  at  grain  level,  and 
texture  development  during  deformation. 

Since  in  crystal  plasticity  models  the  distribution  of  crystal  orientations,  the  available  slip 
and/or  twinning  deformation  systems,  and  the  stress  levels  necessary  to  activate  them  are 
taken  into  account  explicitly,  the  evolution  of  anisotropy  due  to  texture  development  can 
be  characterized  by  measuring  the  initial  texture  and  calculating  grain  reorientation  (i.e. 
updating  the  texture)  using  a  suitable  homogenization  scheme.  Recently,  the  development 
of  crystal  plasticity  models  for  hep  metals  and  the  incorporation  of  crystal  plasticity 
calculations  directly  into  finite  element  (FE)  codes  have  received  much  attention.  Models 
that  account  for  both  slip  and  twinning  activity  and  employ  Taylor  (Bishop  and  Hill, 
1951)  or  self-consistent  averaging  schemes  (Lebensohn  and  Tome,  1993)  to  predict  the 
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aggregate  behavior  have  been  proposed  (e.g.  Staroselsky  and  Anand,  2003;  Salem  et  al., 
2005;  Tome  et  al.  2001,  etc.).  So  far,  these  models  have  been  validated  only  for  quasi¬ 
static  deformation  processes. 

This  paper  presents  a  macroscopic  elastic-viscoplastic  model  for  describing  the  dynamic 
anisotropic  plastic  response  of  textured  metals.  The  aim  is  to  develop  a  model  that  can  be 
applicable  to  both  low-symmetry  and  high-symmetry  pre-textured  materials.  Key  in  this 
development  is  the  use  of  a  recently  proposed  anisotropic  yield  function  (Cazacu  et  al., 
2006).  This  yield  function  (denoted  in  what  follows  as  CPB06)  is  capable  of  describing 
simultaneously  anisotropy  and  tension/compression  asymmetry.  Moreover,  if  the 
parameter  associated  with  strength-differential  effects  is  set  to  zero,  the  CPB06  criterion 
also  describes  accurately  the  yielding  of  metals  with  no  tension/compression  asymmetry. 

The  plan  of  the  paper  is  as  follows.  Section  2  presents  an  overview  of  the  proposed 
elastic-viscoplastic  model.  In  section  3,  numerical  aspects  related  to  the  FE 
implementation  of  the  model  are  discussed.  In  subsections  4.1  and  4.2  we  present  an 
application  of  the  proposed  formulation  to  hep  materials.  In  particular,  we  describe  how 
the  parameters  involved  in  the  evolution  laws  for  a  high-purity  zirconium  clock-rolled 
plate  can  be  determined,  based  on  experimental  information  (texture,  quasi-static  uniaxial 
tension  and  compression  loading  curves  and  metallographic  evidence)  and  simulation 
results  obtained  with  the  viscoplastic  self-consistent  (VPSC)  polycrystalline  model  of 
Lebensohn  and  Tome  (1993).  Next,  a  validation  of  the  anisotropic  low  strain-rate 
formulation  is  provided  by  comparing  FE  simulations  with  experimental  results  on  the 
zirconium  plate.  The  elastic-viscoplastic  extension  of  the  model  is  then  used  to  simulate 
Taylor  cylinder  impact  test  results  for  the  same  material.  We  conclude  with  an 
application  of  the  model  to  bcc  materials  performing  the  analysis  of  the  Taylor  cylinder 
impact  test  on  pre-textured  tantalum.  It  is  shown  that  the  model  describes  with  accuracy 
the  striking  difference  between  the  geometries  of  the  post-test  specimens  of  zirconium 
and  tantalum. 
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2.  ANISOTROPIC  ELASTIC-VISCOPLASTIC  MODEL 


Our  objective  is  to  develop  a  macroscopic  anisotropic  elastic-viscoplastic  model  that 
describes  simultaneously  the  influence  of  strain-rate,  temperature,  and  evolving  texture 
on  the  inelastic  response  of  textured  metals  for  monotonic  loading  paths.  The  basic 
assumption  is  that  the  viscous  properties  of  materials  become  manifest  only  after  the 
passage  to  the  plastic  state.  Thus,  the  strain-rate  £  can  be  decomposed  additively  into  an 
elastic  part  (ee)  and  a  viscoplastic  part  (£vp),  the  latter  representing  the  combined 

viscous  and  plastic  effects,  i.e. 

£  =  «E+£Vp  (1) 

The  evolution  of  the  viscoplastic  strain-rate  is  considered  to  be  given  by  an  overstress 
type  law  of  the  form: 

£»p  =  y(<p(0)|^  (2) 


where 


for  f  <  0, 
for  f  >  0. 


(3) 


In  Eq.  (2),  /  is  the  quasi-static  yield  function,  g  is  the  quasi-static  plastic  potential,  y  a 
viscosity  parameter,  and  m  is  a  strain-rate  sensitivity  constant.  The  quasi-static  yield 
function  is  the  CPB06  yield  function  (Cazacu  et  al.,  2006): 


f^Eyp) 


P(^Evp)  x 
Y(VT) 


(4) 


where  Y  is  a  hardening  parameter,  which  is  a  function  of  the  accumulated  equivalent 
viscoplastic  strain  e  ,  and,  possibly,  of  the  temperature  T,  while  a  is  an  equivalent 

stress  defined  in  the  CPB06  framework  as  (Plunkett  et  al,  2006): 


with  B  being  a  constant  defined  such  that  o  reduces  to  the  tensile  yield  stress  in  the 
rolling  direction,  i.e. 


r+(i*2i-**2r 


and 


a 

(6) 


In  Eq.  (5),  Z,,Z2,S3  are  the  principal  values  of  the  tensor  S  =  L :  S  (where  L  is  a  fourth- 

range  orthotropic  tensor  that  reflects  the  plastic  anisotropy  of  the  material  and  S  is  the 
deviator  of  the  Cauchy  stress  tensor),  k  is  a  strength-differential  parameter,  and  a  is  a 
homogeneity  constant.  Note  that  for  the  case  of  a  =  2,  k  =  0,  and  L  equal  to  the  fourth- 

range  identity  tensor,  CPB06  reduces  to  the  Von  Mises  criterion  and  B  reduces  to  ^3/2  , 

i.e.  the  constant  associated  with  the  von  Mises  effective  stress.  The  effective  plastic  strain 
associated  with  this  anisotropic  yield  function  is  calculated  using  the  principle  of 
equivalence  of  plastic  work  (Hill,  1987). 

The  evolution  of  anisotropy  due  to  evolving  texture  is  accounted  for  by  means  of  a  multi¬ 
scale  procedure,  which  essentially  consists  in  using  experimental  textures  and  stress- 
strain  curves,  polycrystalline  calculations,  and  macroscopic  scale  interpolation 
techniques.  The  procedure  starts  by  fitting  the  parameters  of  a  polycrystal  model  (we  use 
the  VPSC  formulation)  to  reproduce  the  mechanical  response  of  the  pre-textured  material 
along  a  given  deformation  path  (note  that  the  measured  initial  texture  is  needed  as  input 
of  the  polycrystal  model).  Next,  using  the  VPSC  code,  the  evolution  of  texture  is 
calculated  for  fixed  deformation  levels  along  the  given  deformation  path.  Then,  in 

order  to  quantify  the  induced  anisotropy,  the  pre-strained  polycrystal  is  further 
“numerically”  probed  along  various  directions.  The  calculated  yield  stresses  along  these 
various  test  directions  are  in  turn  used  to  determine  the  coefficients 
L(e^),  kCejp),  a(e^)  involved  in  the  macroscopic  CPB06  yield  criterion  and  thus  to 


6 


obtain  the  macroscopic  yield  surface  corresponding  to  the  texture  for  the  given  level  of 
pre-strain,  by  calculating  aJ  according  to  Eq.(5),  as  well  as  YJ  =  Y(sjp).  This  procedure 

is  repeated  for  a  finite  set  of  pre-strain  levels  e^<  8^<...<e^,  j  =  l...n.  Further,  an 

interpolation  procedure  is  used  to  obtain  the  macroscopic  yield  surfaces  corresponding  to 
any  pre-strain  level  (Plunkett  et  al.,  2006),  i.e. 


(8a) 


Y(evp,T)  =  ^(evp)- Yj  +  (!-«£,))•  Y*1 


(8b) 


for  any  <  e^1 ,  j  =  l...m-l.  The  interpolation  parameter  ^(e^) appearing  in  Eqs. 

(8)  is  defined  as: 


—  i+l 


(9) 


It  is  worth  noting  that,  although  the  above  procedure  delivers  a  viscoplastic  response  for 
any  applied  state  at  any  accumulated  deformation,  such  mechanical  response  represents  a 
good  approximation  of  the  actual  behavior  of  the  material  only  for  applied  loadings  with 
moderate  deviations  with  respect  to  the  adopted  pre-strain  deformation  path. 

3.  INTEGRATION  ALGORITHM 

In  order  to  simulate  high-rate  deformation  processes,  the  elastic-viscoplastic  model 
described  by  Eqs.  (1-8)  was  implemented  in  the  commercial  FE  code  ABAQUS  (2003). 
The  algorithmic  aspects  related  to  the  FE  implementation  are  presented  in  what  follows. 

In  displacement-based  FE  formulations,  stress  updates  take  place  at  the  Gauss  points  for  a 
prescribed  nodal  displacement.  We  start  from  time  t,  with  the  known  converged  state 

and  calculate  the  corresponding  values  at  time  t  +  At: 
[ e " Al  ,£typAl  ,o ,+Al ,  s^At  ] .  In  this  incremental  process,  the  total  strain  increment  Ae  is 
decomposed  into  an  elastic  Aee  part  and  a  viscoplastic  part  A£vp  according  to: 
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(10) 


Ae  =  Aee+Aevp 

The  stress  increment  is  related  to  the  elastic  strain  by  Hooke’s  law 
Ao  =  C:AEE=C:(A£-A£vp)  (11) 

where  C  is  the  fourth-range  stiffness  tensor.  Next,  a  trial  elastic  stress:  =  <r'  +  C :  Ae 


is  calculated.  If  f(otria|,evp)^0,  the  increment  is  elastic  and  ot+AI  =  otriiil .  If 

f(o(ri,,,£Vp)>0,  then  there  is  viscoplastic  flow  and  the  stress  update  is  obtained  by 

estimating  the  viscoplastic  strain  increment  A£vp  (Simo  and  Hughes,  1998).  A  fully 

implicit  integration  scheme  was  used,  i.e.  the  viscoplastic  flow  was  determined  by  the 
gradient  of  the  potential  function  at  t  +  At : 


AevpM7A.  =  Y(<r‘)A,(§L 


(12) 


Next,  Eq.  (12)  is  cast  into  the  following  form  (Wang  et  al.,  1997): 


Ac 


vp 


-Ahf® 


^  da  j 


(13) 


t  + At 


where  the  multiplier  Ah  was  introduced  under  the  constraint  that  the  residual  Ar  given 
by: 


A,^.^,  h-) 


(14) 


tends  to  zero  upon  convergence  of  a  local  iterative  procedure.  Since  the  effective  stress  a 
is  homogeneous  of  degree  one  in  stresses,  Aevp=  Ah  .  Substituting  Eq.  (13)  into  Eq.  (11) 

yields: 


Ac  = 


( 


C: 


Ae  -  Ah 


t  +  At  > 


(15) 


To  compute  Ah  ,  a  Taylor  expansion  of  the  residual  is  considered  around  its  value  at  the 
previous  local  iteration,  i.e.: 


8 


(16) 


C“*C"  + 


'  dr  ' 

t+At 

dr 

t+At 

Sr 

8h_+i  + 

:  8<»  , .  + 

_5Ah_ 

n+l 

n 

da 

n+l 

n 

_8T_ 

-it+At 


ST„,«0 


where  n  is  a  counter  for  the  local  iteration  (n  =  0  denotes  the  trial  elastic  state), 
8hn+1=Ahn+,-Ahn,  5<Tn+l=A<Tn+rA<rn,  and  8Tn+1=  ATn+1- ATn.  Replacing  Eq.  (14)  in 
Eq.  (16)  we  obtain  (omitting  supra-indices): 


rn+1  »  rn  +— —  8hn+1  - 

n+l  n  A  .  n+l 

yAt 


S<p 

8a  — 

0tp 

8h  +1- 

0tp 

_da  _ 

n+l 

n 

_5h. 

n+l 

n 

ST 

ST„.,»0 


(17) 


The  stress  update  for  the  trial  state  (n  =  0)  is  found  by  assuming  that  the  total  strain 
increment  was  elastic,  while  for  n  >  0,  the  variation  of  \a  can  be  determined  by 
considering  the  effect  of  the  viscoplastic  strain  on  the  state  of  stress  as: 
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where  ©  ( 0  <  0  <  1 )  is  an  interpolation  parameter.  Combining  Eqs.  (18-19): 
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Further  substitution  of  Eq.  (20)  into  Eq.  (17)  yields 
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In  Eq.  (22)  note  that,  since  the  total  strain  increment  is  fixed  during  the  local  iterative 
procedure,  6£n+1  =  0.  The  derivatives  appearing  in  the  above  equations  are  calculated 
based  on  the  flow  rule  (see  Eq.  2): 
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4.  RESULTS 

4.1  Quasi-static  deformation  of  zirconium 

In  this  subsection  and  the  following  we  present  applications  of  the  proposed  model  to 
high-purity  zirconium.  Results  of  uniaxial  compression  tests  on  a  clock-rolled  zirconium 
plate  were  reported  by  Maudlin  et  al.  (1999b),  Kaschner  and  Gray  (2000),  and  Tome  et 
al.  (2001).  This  hep  aggregate  is  highly  anisotropic,  both  at  single  crystal  and  polycrystal 
levels.  It  was  processed  through  a  series  of  clock-rolling  and  annealing  cycles  to  produce 
a  plate  with  strong  basal  texture  (<c>-axes  of  the  hep  crystals  predominantly  oriented 
along  the  plate  normal  direction).  The  process  of  multiple  rolling  passes  with  rotation  was 
used  in  order  to  obtain  a  nearly  isotropic  in-plane  texture.  Let  (x,  y,  z)  be  a  coordinate 
reference  frame  associated  with  the  plate,  where  the  (jc,  y)  plane  corresponds  to  the 
original  rolling  plane,  while  the  z-direction  is  the  through-thickness  (TT)  direction.  Right- 
circular  cylindrical  tests  specimens  were  cut  from  the  plate  both  in  the  TT  (z-)  and  in¬ 
plane  (IP)  ( x -  or  y-)  directions.  Quasi-static  tests  in  both  IP  and  TT  directions  were 


10 


carried  out  at  room  temperature  (298  K)  and  liquid  nitrogen  temperature  (77  K)  (Tome  et 
al.,  2001).  In  this  paper  we  will  restrict  our  analysis  to  room  temperature  results.  It  was 
found  that  the  mechanical  response  is  highly  dependent  on  the  predominant  orientation  of 
the  <c>-axes  with  respect  to  the  loading  direction  (Tome  et  al.,  2001).  The  large 
variations  observed  in  flow  stress  and  strain-hardening  rate  are  due  to  the  relative 
orientation  of  the  available  twinning  and  slip  systems  and  the  critical  stress  levels 
required  to  activate  them.  In  the  case  of  TT  samples,  compressive  loading  was  applied 
perpendicular  to  the  basal  plane,  thus  the  hard  pyramidal  <c+a>-slip  was  the  dominant 
mechanism.  Compressive  loading  of  the  IP  samples  was  favorable  to  activate  prismatic 
slip  and  tensile  twinning,  therefore  plastic  flow  occurred  at  lower  applied  stresses.  In 
what  follows,  all  the  reported  numerical  tests  using  the  VPSC  model  were  carried  out 
assuming:  a)  the  experimental  initial  texture  (consisting  of  377  orientations)  of  the  clock- 
rolled  zirconium  plate,  b)  the  deformation  mechanisms  characterized  as  operational  at 
room  temperature  (i.e.  prismatic  <a>-slip,  pyramidal  <c+a>-slip  and  tensile  twinning), 
and  c)  the  values  of  the  slip  and  twinning  parameters  (critical  stresses,  hardening 
coefficients,  rate  sensitivity  exponent),  as  reported  by  Tome  et  al.  (2001). 

As  already  mentioned,  in  order  to  experimentally  characterize  the  anisotropy  induced  by 
texture  evolution,  the  material  should  be  subjected  to  a  specific  loading,  up  to  pre¬ 
determined  deformation  levels.  The  choice  of  the  pre-strain  loading  direction  should  be 
consistent  with  the  main  deformation  path  that  the  aggregate  will  undergo  in  a  specific 
application  (e.g.  uniaxial  compression  for  Taylor  impact  testing).  Then,  specimens  should 
be  cut  from  the  pre-strained  sheet  and  further  tests  performed  to  quantify  the  anisotropy 
induced  by  deformation.  However,  this  would  require  a  large  number  of  mechanical  tests, 
which  could  be  very  difficult  to  carry  out  because,  e.g.  at  large  compressive  pre-straining 
levels  buckling  of  the  sheet  is  likely  to  occur.  Alternatively,  the  evolution  of  the  yield 
surface  can  be  studied  by  performing  numerical  tests  using  VPSC  model.  For  the  clock- 
rolled  plate  of  zirconium  described  previously,  the  VPSC  code  was  first  used  to 
characterize  the  evolution  of  the  yield  surface  during  IP  compression.  For  this  purpose, 
the  polycrystal  was  pre-strained  up  to  given  deformation  levels.  Figure  1  shows  the 
biaxial  stress  plane  projections  of  the  VPSC  yield  loci  (symbols)  corresponding  to  pre- 
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strain  levels  of  =0.2%,  e^=l%,  e^  =  5%,  e^,  =25%,  =35%,  =45%,  and 

ev7p  =60%.  Next,  for  each  individual  pre-strain  level  we  obtained  oj  and  determined  the 

coefficients  L(e^,),  k(e^),  a(e^)  (see  Eqs.  5-7),  and  Yj=Y(eJp)  as  well.  The 

homogeneity  parameter  was  considered  to  have  a  fixed  value  a  =  2  because  of  the  nearly 
elliptical  shape  of  the  yield  loci  for  high-purity  zirconium.  This  procedure  was  repeated 
for  all  seven  pre-strain  levels.  The  numerical  values  of  these  coefficients  for  each  CPB06 
yield  surface  are  listed  in  Table  1.  These  coefficients  quantify  the  very  strong  plastic 
anisotropy  for  this  material  and  its  strong  evolution  with  strain.  Also,  the  linear 
interpolation  procedure  was  used  to  obtain  the  macroscopic  yield  surfaces  corresponding 
to  any  pre-strain  level  (see  Eq.  8).  In  Fig.l,  the  dashed  lines  represent  the  CPB06  yield 
surfaces  for  intermediate  pre-strain  levels  obtained  using  linear  interpolation.  Note  that  an 
accurate  description  of  the  evolving  yield  function  for  pre-strain  levels  of  around  30% 
and,  further,  up  to  60%  are  needed  for  the  quasi-static  IP  compression  simulation  below, 
and  for  the  Taylor  cylinder  simulation  of  next  subsection,  respectively.  It  is  also  worth 
mentioning  that  the  successive  yield  function  projections  for  pre-strains  higher  than  5% 
are  more  rounded  and  centered,  reflecting  a  decrease  of  the  plastic  anisotropy  (note  that 
all  anisotropy  coefficients  have  initially  rather  large  values  which  decrease  at  large 
strains).  This  effect  is  related  to  the  high  tensile  twinning  activity  and  the  concomitant 
strong  reorientation  of  the  twinned  portions  of  the  grains  that  take  place  during  in-plane 
compression. 

The  proposed  model  was  used  to  simulate  the  uniaxial  compression  of  a  right  cylinder  of 
circular  cross  section  with  its  axis  parallel  to  the  in-plane  jc-direction  up  to  28% 
longitudinal  strain,  for  comparison  with  analog  experimental  results  reported  by  Tome  et 
al.  (2001)  (Fig.  2).  Note  that  because  the  z-direction  is  very  hard  to  deform,  the  specimen 
deforms  mostly  in  the  _y-direction  and  the  cross-section  of  the  deformed  specimen 
becomes  elliptical.  Our  simulations  predicts  in-plane  expansion  of  22%  and  out-of-plane 
expansion  of  6%  (or  an  ovalization  long/short  radii=1.19),  which  agrees  very  well  with 
the  experimental  results  reported  in  Tome  et  al.  (2001).  It  is  worth  noting  that  under  the 
assumption  of  isotropic  hardening  and  for  fixed  values  of  the  anisotropy  coefficients  (i.e. 
texture  evolution  was  neglected  and  the  anisotropy  coefficients  for  any  pre-strain  level 
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were  set  equal  to  the  values  of  L  obtained  based  on  VPSC  yield  stresses  calculated  by 
probing  the  0.2%  pre-strained  polycrystal  in  different  directions),  the  predicted 
expansions  would  have  been  of  27%  in-plane  and  1%  out-of-plane.  This  clearly  implies 
that  a  good  agreement  with  the  experimental  results  can  be  obtained  only  if  the  changes 
in  anisotropy  induced  by  texture  evolution  (in  this  case,  a  decrease  in  plastic  anisotropy) 
is  taken  into  account.  Figure  2a  shows  the  final  ABAQUS  FE  mesh.  The  photograph  of 
the  final  experimental  cross-section  shape  and  the  superimposed  simulation  result 
(obtained  using  the  FE  program  EPIC  directly  linked  to  the  VPSC  code)  are  shown  in 
Fig.  2b  (simulation  and  data  after  Tome  et  al,  2001).  The  comparison  demonstrates  that 
very  good  agreement  with  the  experimental  results  can  be  obtained  by  using  an 
appropriate  initial  yield  surface  representation  (CPB06)  and  a  hardening  law  that 
accounts  for  texture  evolution. 

Next,  we  have  simulated  the  quasi-static  through-thickness  compression  of  the  zirconium 
specimen.  Figure  3  shows  the  yield  stresses  obtained  with  the  VPSC  code  (symbols)  for 
five  individual  levels  of  compressive  pre-strain  (0.2%,  1%,  5%,  25%  and  35%)  along  the 
through-thickness  direction,  the  corresponding  five  individual  CPB06  yield  surfaces  for 
which  the  coefficients  were  determined  (solid  lines),  and  yield  surfaces  at  intermediate 
levels  of  strain  obtained  by  linear  interpolation  (dashed  lines).  Note  that  a  description  of 
the  evolving  yield  function  up  to  pre-strain  levels  of  around  35%  was  enough  for  the  28% 
strain  quasi-static  TT  compression  simulation  below. 

The  numerical  values  of  the  anisotropy  coefficients  for  CPB06  corresponding  to  each 
individual  level  of  pre-strain  are  listed  in  Table  2.  Note  that  the  CPB06  yield  surface 
evolution  associated  to  deformation  induced  anisotropy  in  uniaxial  compression  is  very 
different  than  the  evolution  associated  to  uniaxial  compression  in  the  in-plane  direction 
(cf.  Figs.  1  and  3).  Through-thickness  compression  is  dominated  by  slip  so  the  texture 
does  not  significantly  change  during  loading,  whereas  during  in-plane  compression 
tensile  twinning  is  easily  activated  thus  producing  a  strong  evolution  in  texture. 

The  yield  surface  evolution  depicted  in  Fig.  3  was  in  turn  used  in  an  ABAQUS  FE 
simulation  of  compression  of  a  right  cylinder  of  circular  cross  section  and  axis  parallel  to 
the  z-axis  up  to  28%  longitudinal  strain,  and  compared  with  the  experimental  results 
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(Tome  et  al,  2001)  (Fig.  4).  Since  compression  takes  place  along  the  axis  of  symmetry  of 
the  sample,  the  final  cross-section  remains  circular.  The  proposed  model  correctly 
predicts  this  experimental  trend  (cf.  Figs.  4a  and  4b). 


4.2  High  strain-rate  deformation  of  zirconium 

Controlled  impact  tests  are  used  to  characterize  the  high  strain-rate  mechanical  response 
of  metals.  Such  tests  provide  results  in  the  1 04- 1 05  s'1  strain-rate  range,  directly  above 
that  data  accessible  via  Hopkinson  Bar  techniques  and  below  shock  experimentation.  One 
standard  approach  used  is  the  impact  of  a  right  circular  cylinder  on  either  a  massive  anvil 
(rigid  plate)  or  a  receptor  rod  of  identical  material.  This  experiment  is  commonly  called  a 
Taylor  impact  test  (Taylor,  1948).  Due  to  the  impact,  the  rod  plastically  deforms  and 
shortens  causing  the  material  at  the  impact  surface  to  flow  radially  outward  relative  to  the 
rod  axis. 

Maudlin  et  al.  (1999b)  reported  results  of  Taylor  cylinder  impact  tests  using  in-plane 
specimens  cut  from  the  same  zirconium  plate  described  and  characterized  in  the  previous 
section  (i.e.  cylindrical  axes  were  coincident  with  in-plane  directions).  The  specimens 
were  50.8  mm  in  length  and  have  a  length  over  diameter  ratio  of  6.67.  The  anvil  was 
AF1410  steel,  heat-treated  to  a  surface  hardness  of  Rc  58  and  lapped  to  a  mirror  finish. 
The  cylinders  were  fired  using  a  gas-driven  gun  at  various  velocities:  50,  101,  170,  and 
243  m/s.  Only  post-test  data  for  the  243  m/s  impact  velocity  were  reported.  As  discussed 
in  Maudlin  et  al.  (1999b),  the  data  set  consists  of  high-resolution  scans  of  the 
photographed  footprint,  i.e.  the  cross-sectional  area  at  the  impact  interface  (on  the  y-z 
plane)  and  of  the  photographed  side  profiles  for  the  major  and  minor  dimensions.  From 
the  analysis  of  the  digitized  side  profiles,  axial  strain  profiles  (log  strain)  as  a  function  of 
axial  position,  measured  relative  to  the  impact  interface,  were  obtained. 

We  have  applied  the  proposed  elastic-viscoplastic  model  to  the  description  of  the  high 
strain-rate  deformation  of  the  clock-rolled  zirconium  plate  at  room  temperature  and 
compared  the  simulations  with  Taylor  impact  test  results  reported  by  Maudlin  et  al 
(1999b).  In  addition  to  the  yield  surface  evolution  due  to  plastic  deformation,  the 
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viscosity  coefficient  y  and  the  strain-rate  sensitivity  parameter  m  (see  Eqs.  2-3)  needed 
to  be  determined.  Results  of  dynamic  in-plane  tests  performed  on  the  clock-rolled  Zr 
material  at  strain-rates  from  1000/s  to  3500/s  conducted  using  a  Split  Hopkinson  Pressure 
Bar  (SHPB)  were  reported  by  Kaschner  and  Gray  (2000).  Based  on  these  tests,  we  have 
determined  y  =  2500  s-1  and  m  =  7.0  (see  Fig.  5).  The  constitutive  behavior  for  the  in¬ 
plane  compression  corresponding  to  strains  beyond  the  range  of  the  experimental  data 
was  estimated  based  on  polycrystalline  calculations. 

In  order  to  simulate  Taylor  impact  tests  of  zirconium,  the  proposed  elastic-viscoplastic 
model  was  implemented  into  an  ABAQUS  explicit  user  material  subroutine.  The  rate- 
independent  yield  surface  characterization  was  the  one  described  in  the  previous 
subsection,  based  on  VPSC  calculations  of  the  induced  anisotropy  during  quasi-static  in¬ 
plane  compression  coupled  with  CPB06  (see  Fig.  3).  It  should  be  noted  that  by  making 
this  choice,  we  have  disregarded  the  possibility  of  changes  in  the  combination  of  active 
deformation  mechanisms  at  single  crystal  level  as  strain-rate  increases.  For  comparison 
purposes,  simulations  under  the  assumption  of  isotropic  hardening  and  anisotropy 
coefficients  set  equal  to  the  anisotropy  coefficients  corresponding  to  0.2%  pre-strain 
(texture  evolution  neglected)  were  also  carried  out. 

The  zirconium  cylinder  was  modeled  using  2117  ABAQUS  C3D8R  linear  brick  elements 
with  free  boundary  conditions.  The  anvil  target  was  modeled  as  an  analytical  rigid 
surface.  Simulations  were  done  for  the  impact  velocity  of  243  m/s.  The  step  size  used  to 
simulate  the  event  was  At  =2-10'8s .  Calculations  were  performed  for  90  ps,  at  which 
time  the  specimen  had  rebound  off  of  the  target  and  all  plastic  deformation  had  ceased. 
Figure  6  shows  a  comparison  between  the  experimental  and  simulated  logarithmic  strain 
profile  along  the  major  and  minor  axes  of  the  post-test  specimen.  The  strains  along  the 
major  and  minor  elliptical  axes,  which  are  both  zero  at  large  z,  accumulate  at  different 
rates  as  z  decreases.  A  good  agreement  with  the  experimental  data  is  observed.  Figure  7 
depicts  a  visual  comparison  between  the  simulated  and  experimental  major  and  minor 
profiles  and  footprint  of  the  post-test  specimen.  The  slight  differences  between  the 
experimental  data  and  the  simulation  results  are  most  likely  related  to  temperature  effects 
which  were  not  taken  into  consideration  in  the  simulations  because  experimental  data  was 
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not  available.  Note  that  the  geometry  of  the  footprint  of  the  specimen  clearly  indicates  a 
decrease  in  anisotropy  (the  ratio  of  major  to  minor  diameter  is  approximately  1.1).  This 
is  due  to  the  large  twinned  volume  fraction  resulting  from  strains  of  around  60-70%  at  the 
impact  section  (see  strain  levels  in  Fig.  7)  which  tends  to  decrease  the  intensity  of  the 
texture  (as  already  seen  in  Fig.l  and  discussed  in  subsection  4.1). 

As  discussed  earlier  in  the  quasi-static  case,  the  consideration  of  a  hardening  law  that 
accounts  for  texture  evolution  is  essential  to  obtain  the  good  agreement  with  the 
experimental  results.  Figure  8  shows  a  comparison  between  experimental  minor  and 
major  profiles  of  the  post-test  specimen  and  those  calculated  using  isotropic  hardening 
and  fixed  values  for  the  anisotropy  coefficients  corresponding  to  0.2%  pre-strain.  It  is 
evident  that  the  effects  of  substantial  twinning  activity  can  be  captured  only  by  a 
hardening  law  that  reflects  the  evolution  of  texture. 


4.3  High  strain-rate  deformation  of  tantalum 

The  anisotropic  elastic-viscoplastic  model  was  originally  developed  to  capture  the 
anisotropic,  asymmetric  hardening  behavior  of  hep  metals.  The  static  yield  surface  in  this 
model  is  the  CPB06  yield  surface.  However,  the  CPB06  criterion  can  be  also  used  to 
represent  the  yielding  behavior  of  a  material  with  cubic  crystallographic  structure 
deforming  solely  by  slip  by  simply  setting  the  strength-differential  coefficient  k  to  zero. 
In  this  section,  we  apply  the  proposed  formulation  to  the  description  of  the  anisotropy  of 
a  bcc  tantalum  sheet  (obtained  by  unidirectional  rolling  and  displaying  a  near-orthotropic 
texture)  subjected  to  high  strain-rate  and  large  plastic  deformations  (data  after  Maudlin  et 
al.,  1 999a, b). 

Unlike  for  the  previously  studied  clock-rolled  zirconium,  in  this  case  there  was  no 
sufficient  experimental  data  available  to  quantify  the  deformation  induced  anisotropy 
using  VPSC  simulations  (i.e.  no  data  based  on  which  we  could  determine  the  parameters 
involved  in  the  VPSC  model  were  available).  Therefore,  to  estimate  the  coefficients 
involved  in  the  macroscopic  CPB06  yield  criterion,  we  have  used  the  values  of  yield 
stresses  in  different  directions  reported  by  Maudlin  et  al.  (1999a,b).  The  numerical  values 
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of  the  CPB06  coefficients  are  listed  in  Table  3.  Moreover,  due  to  lack  of  information 
regarding  the  yield  behavior  in  pure  shear,  it  was  assumed  a  value  equal  to  0.577  times 
the  yield  stress  in  uniaxial  tension  along  the  rolling  direction.  Figure  9  shows  the 
projection  in  the  biaxial  plane  of  the  CPB06  orthotropic  yield  surface  for  tantalum  (solid 
lines)  together  with  the  yield  stresses  (symbols)  calculated  by  Maudlin  et  al.  (1999a,b) 
using  a  Taylor-Bishop-Hill  (TBH)  polycrystal  model. 


Uniaxial  compression  test  results  at  various  strain-rates  and  temperatures  were  reported 
for  tantalum  specimens  cut  from  the  same  rolled  plate.  Results  of  Split-Hopkinson  bar 
tests  performed  at  different  temperatures  in  the  range  200°C-1000°C  and  for  strain-rates 
in  the  range  2800-3900  s'1  have  shown  that  the  compressive  response  depends  on  the 
strain-rate  and  temperature.  The  yield  stress  was  found  to  increase  with  increasing  strain- 
rate  and  decreasing  temperature.  To  capture  the  effect  of  temperature  on  the  compressive 
yield  stress  of  tantalum  we  used  a  Johnson-Cook  type  hardening  law  (Johnson  and  Cook, 
1983): 


Y(vT)  =  Y.<S) 


T-T. 


x  tm 


T  -T 

V  1melt  Arm  J 


(24) 


In  Eq.  (24),  Y^e^)  is  the  rate-independent  hardening  law  at  room  temperature,  Tm,  is  the 

room  temperature  (typically  25°C),  Tmc|t  is  the  melting  temperature,  and  tm  is  a  material 
constant.  For  tantalum,  Tmelt  =  2977°C,  while  a  value  of  0.42  for  the  exponent  tm  gives 

an  overall  good  approximation  of  the  experimental  data.  Assuming  that  all  mechanical 
work  is  converted  into  heat,  from  the  energy  equation  (Malvern,  1969)  the  temperature 
rise  due  to  mechanical  deformation  is: 


a :  \svn 
AT  = - ^ 

Pc„ 


(25) 


For  tantalum,  p=  16640  kg/m3 ,  cp=140  J/kgK  .  Furthermore,  using  the  available  data,  a 

viscosity  /  =  1500s'1  and  a  strain-rate  sensitivity  parameter  m  =  5.0  (see  Eq.  2)  were 
determined. 
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Next,  we  evaluated  the  predictive  capability  of  the  present  model  to  simulate  the  severe 
geometry  changes  of  the  tantalum  Taylor  specimens  reported  by  Maudlin  et  al.  (1999a,b). 
The  cylinders  were  38.1  mm  long  and  7.6  mm  in  diameter.  The  longitudinal  axis  of  the 
specimen  was  either  parallel  to  the  rolling  direction  (y )  or  transverse  direction  (x).  Only 
results  on  specimens  having  the  longitudinal  axis  parallel  to  the  rolling  direction  were 
reported.  The  cylinders  were  launched  using  a  powder  gun.  The  velocity  of  the  cylinders 
measured  using  both  pressure  transducers  and  parallel  laser  beams  were  1 75  m/s  ±  3  m/s. 
The  anvil  target  was  made  of  steel.  Metallographic  analysis  of  the  pre-  and  post-test 
specimens  showed  negligible  texture  evolution.  Therefore,  the  experiment  was  modeled 
assuming  isotropic  hardening.  Because  the  material  was  orthotropic,  only  a  quarter 
section  was  modeled  using  8170  ABAQUS  C3D8R  linear  brick  elements  with  free 
boundary  conditions  and  an  initial  velocity  of  175  m/s.  The  step  size  used  to  simulate  the 
event  was  At=2xl0'8s  and  was  conducted  until  100  ps.  At  100  ps  the  specimen  had 
rebound  off  of  the  target,  and  all  plastic  deformation  had  ceased.  The  anvil  target  was 
modeled  as  an  analytical  rigid  surface.  A  higher  number  of  elements  were  used  for  the 
tantalum  simulations  as  compared  to  the  zirconium  simulations  because  the  tantalum 
specimen  undergoes  larger  deformations  (for  tantalum  at  z  =  0,  the  strain  along  the  major 
and  minor  dimensions  were  of  87%  and  70%,  respectively  while  for  zirconium  the  values 
of  the  strains  at  z=0  were  of  44%  and  31%,  respectively).  Axial  strain  profiles  for  both 
the  major  and  minor  sides  of  the  deformed  specimens  (cylinder  axis  parallel  to  the  rolling 
direction)  in  comparison  with  experiments  are  shown  in  Fig.  10.  Very  good  agreement  is 
observed.  The  simulated  and  experimental  footprints  of  the  post  test  specimen  are  shown 
in  Fig.  11.  The  data  shown  in  Figs.  10-11  are  the  results  from  three  separate  experiments 
conducted  for  different  specimens  under  the  same  conditions.  The  simulated  and 
experimental  post-test  specimens  are  visually  compared  in  Fig.  12.  In  Maudlin  et  al. 
(1999b),  the  same  tests  were  simulated  by  introducing  anisotropy  through  direct  link 
between  the  Taylor-Bishop-Hill  polycrystal  model  and  the  FE  code  EPIC,  assuming  the 
presence  of  a  polycrystal  at  each  integration  point.  Their  results  were  also  in  very  good 
agreement  with  the  experimental  data.  However,  Figs.  10-12  clearly  demonstrate  that 
very  accurate  results  can  also  be  obtained  with  the  present  formulation,  at  a  fraction  of 
the  computation  time. 
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5.  CONCLUSIONS 


A  macroscopic  anisotropic  elastic-viscoplastic  model  that  captures  the  influence  of 
evolving  texture  on  the  mechanical  response  of  textured  metals  for  both  quasi-static  and 
dynamic  loading  conditions  was  proposed.  Yielding  was  described  using  the  anisotropic 
yield  function  proposed  by  Cazacu  et  al  (2006).  This  yield  function  can  capture 
simultaneously  both  anisotropy  and  tension/compression  asymmetry.  Moreover,  if  the 
parameter  associated  with  strength  differential  effects  is  set  to  zero,  the  yield  function 
also  describes  with  accuracy  yielding  of  metals  with  cubic  symmetry. 

The  anisotropy  coefficients  as  well  as  the  size  of  the  elastic  domain  were  considered  to  be 
functions  of  the  accumulated  plastic  strain.  To  capture  the  anisotropy  induced  by  texture 
evolution,  a  multi-scale  methodology  consisting  in  using  experimental  textures  and 
stress-strain  curves,  VPSC  polycrystalline  calculations,  and  macroscopic  scale 
interpolation  techniques  was  used.  The  procedure  started  by  fitting  the  parameters  of  the 
polycrystal  model  to  reproduce  the  mechanical  response  of  the  pre-textured  material. 
Next,  using  this  adjusted  polycrystal  model,  the  evolution  of  texture  was  calculated  for 
fixed  pre-strain  level  along  a  given  deformation  path.  Then,  in  order  to  quantify  the 
induced  anisotropy,  the  pre-strained  polycrystal  was  numerically  probed  along  various 
directions.  The  calculated  yield  stresses  along  these  directions  were  in  turn  used  to 
determine  the  CPB06  coefficients.  This  procedure  was  repeated  for  a  finite  set  of  pre¬ 
strain  levels  and  linear  interpolation  was  in  turn  used  to  obtain  the  macroscopic  yield 
surfaces  corresponding  to  an  arbitrary  pre-strain  level 

The  above  procedure  was  used  to  model  the  quasi-static  mechanical  response  of  high- 
purity  zirconium.  Since  the  proposed  model  accounts  for  the  deformation  mechanisms 
associated  with  a  given  deformation  path,  it  captured  the  strong  anisotropy  and  texture 
evolution  produced  by  twinning  which  is  observed  when  the  material  is  subjected  to  in¬ 
plane  compression.  On  the  other  hand,  the  resulting  macroscopic  yield  surface  evolution 
along  through-thickness  compression  was  very  different  than  the  one  associated  with  in¬ 
plane  compression  (cf.  Figs.  1  and  3).  Through-thickness  compression  is  dominated  by 
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slip  so  the  texture  does  not  significantly  change  during  loading,  whereas  during  in-plane 
compression  tensile  twinning  is  easily  activated,  producing  a  strong  texture  evolution. 

Based  on  the  above  model  for  quasi-static  behavior,  an  overstress  type  of  approach  was 
in  turn  used  to  account  for  rate  effects.  The  resulting  elastic-viscoplastic  model  was 
applied  to  the  description  of  the  high  strain-rate  response  of  low-symmetry  (clock-rolled 
hep  zirconium)  and  high-symmetry  (rolled  bee  tantalum)  textured  metals.  The  very  good 
agreement  between  the  simulated  and  experimental  post-test  geometries  of  the  Taylor- 
impact  specimens  in  terms  of  major  and  minor  side  profiles  and  impact  footprints  proved 
the  ability  of  the  model  to  describe  the  striking  difference  in  the  evolution  of  anisotropy 
between  zirconium  and  tantalum. 
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FIGURE  CAPTIONS 


Figure  1:  Theoretical  yield  surface  evolution  for  a  zirconium  clock-rolled  plate  during 
in-plane  compression  for  various  levels  of  pre-strain.  Solid  lines:  CPB06  yield  surfaces 
corresponding  to  fixed  pre-strain  levels,  determined  using  the  VPSC  code  (symbols). 
Dashed  lines:  CPB06  yield  surfaces  for  intermediate  pre-strain  levels  obtained  by  linear 
interpolation. 

Figure  2:  Cross-sections  for  quasi-static  compression  of  a  zirconium  cylinder  with  axis 
parallel  to  an  in-plane  direction,  (a)  Final  ABAQUS  FE  mesh  obtained  with  the  proposed 
model,  (b)  Photograph  of  the  final  experimental  cross-section  and  EPIC/VPSC  simulation 
results  (dashed  lines),  after  Tome  et  al.  (2001).  The  arrows  indicate  the  predominant 
direction  of  the  crystallographic  <c>-axes. 

Figure  3:  Theoretical  yield  surface  evolution  for  a  zirconium  clock-rolled  plate  during 
trough-thickness  compression  for  various  levels  of  pre-strain.  Solid  lines:  CPB06  yield 
surfaces  corresponding  to  fixed  pre-strain  levels,  determined  using  the  VPSC  code 
(symbols).  Dashed  lines:  CPB06  yield  surfaces  for  intermediate  pre-strain  levels  obtained 
by  linear  interpolation. 

Figure  4:  Cross-sections  for  quasi-static  compression  of  a  zirconium  cylinder  with  axis 
parallel  to  the  through-thickness  direction,  (a)  Final  ABAQUS  FE  mesh  obtained  with  the 
proposed  model,  (b)  Photograph  of  the  final  experimental  cross-section  and  EPIC/VPSC 
simulation  results  (dashed  lines),  after  Tome  et  al.  (2001). 

Figure  5:  Stress-strain  curves  of  in-plane  compression  of  clock-rolled  zirconium  samples 
for  various  strain-rates.  Symbols:  experimental  data,  after  Kaschner  and  Gray  (2000). 
Lines:  simulation  results  using  the  proposed  elastic-viscoplastic  model  for  y  =  2500  s' 
and  m  =  7.0 . 
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Figure  6:  Comparison  of  theoretical  (lines)  and  experimental  (symbols)  logarithmic 
strain  profiles  of  the  post-test  zirconium  Taylor  impact  specimen.  Data  after  Maudlin  et 
al.  (1999b). 

Figure  7:  Comparison  of  the  simulated  (with  contours  of  axial  plastic  strain)  and 
experimental  cross-sections  of  the  post-test  zirconium  Taylor  impact  experiment  for  the 
major  profile,  the  minor  profile,  and  the  footprint.  Data  after  Maudlin  et  al.  (1999b). 

Figure  8:  Comparison  of  theoretical  (lines)  and  experimental  (symbols)  logarithmic 
strain  profiles  of  the  post-test  zirconium  Taylor  impact  specimen.  Simulation  were 
performed  with  an  isotropic  hardening  law.  Data  after  Maudlin  et  al.  (1999b). 

Figure  9:  Line:  CPB06  yield  surface  corresponding  to  a  rolled  sheet  of  tantalum. 
Symbols:  yield  stress  points  predicted  with  the  Taylor-Bishop-Hill  approach.  Stresses  are 
normalized  about  the  yield  stress  in  the  x-direction  (rolling  direction).  Data  after  Maudlin 
et.  al.  (1999a,b). 

Figure  10:  (a)  Major  and  (b)  minor  profiles  of  the  tantalum  Taylor  impact  specimen. 
Solid  lines  represent  simulation  results,  symbols  represent  experimental  data.  Data  after 
Maudlin  et.  al.  (1999a,b). 

Figure  11:  Footprint  of  the  tantalum  Taylor  impact  specimen.  Solid  lines  represent 
simulation  results,  symbols  represent  experimental  data.  Data  after  Maudlin  et.  al. 
( 1 999a,b). 

Figure  12:  Comparison  between  simulated  (with  contours  of  axial  plastic  strain)  and 
experimental  tantalum  Taylor  impact  specimens.  Photographs  of  the  post-test  specimen 
after  Maudlin  et  al.  (1999a,b). 
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TABLE  CAPTIONS 


Table  1:  CPB06  yield  function  coefficients  corresponding  to  various  pre-strain  levels  for 
uniaxial  compression  of  clock-rolled  zirconium  along  an  in-plane  direction.  For  all  levels 
of  effective  plastic  strain:  a  =  2  and  Li  i  =  1 .0 


Table  2:  CPB06  yield  function  coefficients  corresponding  to  various  pre-strain  levels  for 
uniaxial  compression  of  clock-rolled  zirconium  along  the  through-thickness  direction.  . 
For  all  levels  of  effective  plastic  strain:  a  =  2  and  Ln  =  1.0 


Table  3:  CPB06  yield  function  coefficients  corresponding  to  rolled  tantalum;  a  =  2  and 

Ln  =  1.0. 
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Figure  1:  Theoretical  yield  surface  evolution  for  a  zirconium  clock-rolled  plate  during 
in-plane  compression  for  various  levels  of  pre-strain.  Solid  lines:  CPB06  yield  surfaces 
corresponding  to  fixed  pre-strain  levels,  determined  using  the  VPSC  code  (symbols). 
Dashed  lines:  CPB06  yield  surfaces  for  intermediate  pre-strain  levels  obtained  by  linear 
interpolation. 
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(a)  (b) 


Figure  2:  Cross-sections  for  quasi-static  compression  of  a  zirconium  cylinder  with  axis 
parallel  to  an  in-plane  direction,  (a)  Final  ABAQUS  FE  mesh  obtained  with  the  proposed 
model,  (b)  Photograph  of  the  final  experimental  cross-section  and  EPIC/VPSC  simulation 
results  (dashed  lines),  after  Tome  et  al.  (2001).  The  arrows  indicate  the  predominant 
direction  of  the  crystallographic  <c>-axes. 
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Figure  3:  Theoretical  yield  surface  evolution  for  a  zirconium  clock-rolled  plate  during 
trough-thickness  compression  for  various  levels  of  pre-strain.  Solid  lines:  CPB06  yield 
surfaces  corresponding  to  fixed  pre-strain  levels,  determined  using  the  VPSC  code 
(symbols).  Dashed  lines:  CPB06  yield  surfaces  for  intermediate  pre-strain  levels  obtained 
by  linear  interpolation. 
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Figure  4:  Cross-sections  for  quasi-static  compression  of  a  zirconium  cylinder  with  axis 
parallel  to  the  through-thickness  direction,  (a)  Final  ABAQUS  FE  mesh  obtained  with  the 
proposed  model,  (b)  Photograph  of  the  final  experimental  cross-section  and  EPIC/VPSC 
simulation  results  (dashed  lines),  after  Tome  et  al.  (2001). 
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Figure  5:  Stress-strain  curves  of  in-plane  compression  of  clock-rolled  zirconium  samples 
for  various  strain-rates.  Symbols:  experimental  data,  after  Kaschner  and  Gray  (2000). 
Lines:  simulation  results  using  the  proposed  elastic-viscoplastic  model  for  y  =  2500  s"1 
and  m  =  7.0 . 
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Figure  6:  Comparison  of  theoretical  (lines)  and  experimental  (symbols)  logarithmic 
strain  profiles  of  the  post-test  zirconium  Taylor  impact  specimen.  Data  after  Maudlin  et 
al.  (1999b). 
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Figure  7:  Comparison  of  the  simulated  (with  contours  of  axial  plastic  strain)  and 
experimental  cross-sections  of  the  post-test  zirconium  Taylor  impact  experiment  for  the 
major  profile,  the  minor  profile,  and  the  footprint.  Data  after  Maudlin  et  al.  (1999b). 
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Figure  8:  Comparison  of  theoretical  (lines)  and  experimental  (symbols)  logarithmic 
strain  profiles  of  the  post-test  zirconium  Taylor  impact  specimen.  Simulation  were 
performed  with  an  isotropic  hardening  law.  Data  after  Maudlin  et  al.  (1999b). 
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Figure  9:  Line:  CPB06  yield  surface  corresponding  to  a  rolled  sheet  of  tantalum. 
Symbols:  yield  stress  points  predicted  with  the  Taylor-Bishop-Hill  approach.  Stresses  are 
normalized  about  the  yield  stress  in  the  x-direction  (rolling  direction).  Data  after  Maudlin 
et.  al.  (1999a,b). 
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(a)  (b) 


Figure  10:  (a)  Major  and  (b)  minor  profiles  of  the  tantalum  Taylor  impact  specimen. 
Solid  lines  represent  simulation  results,  symbols  represent  experimental  data.  Data  after 
Maudlin  et.  al.  (1999a,b). 
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Figure  11:  Footprint  of  the  tantalum  Taylor  impact  specimen.  Solid  lines  represent 
simulation  results,  symbols  represent  experimental  data.  Data  after  Maudlin  et.  al. 
(1999a,b). 
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Figure  12:  Comparison  between  simulated  (with  contours  of  axial  plastic  strain)  and 
experimental  tantalum  Taylor  impact  specimens.  Photographs  of  the  post-test  specimen 
after  Maudlin  et  al.  (1999a,  b). 
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Table  1:  CPB06  yield  function  coefficients  corresponding  to  various  pre-strain  levels  for 
uniaxial  compression  of  clock-rolled  zirconium  along  an  in-plane  direction.  For  all  levels 
of  effective  plastic  strain:  a  =  2  and  Li  i  =  1.0 


k 

L,2 

L,3 

L22 

L23 

L33 

L44 

L55 

1-66 

0.2 

0.071 

2.883 

1.875 

1.252 

1.965 

1.204 

0.989 

1.075 

3.956 

1% 

0.224 

2.15 

1.407 

1.179 

1.438 

0.892 

0.604 

0.617 

1.691 

5% 

0.500 

2.923 

1.823 

0.865 

1.745 

0.418 

0.984 

1.145 

2.980 

25 

0.263 

2.704 

1.727 

1.821 

2.056 

0.869 

0.949 

0.791 

1.403 

35 

0.174 

2.505 

1.683 

2.320 

2.054 

1.012 

0.889 

0.694 

1.098 

45 

-0.269 

0.853 

1.444 

1.143 

0.624 

1.007 

0.638 

0.461 

0.665 

60 

-0.220 

0.766 

0.858 

0.637 

0.900 

0.644 

0.299 

0.189 

0.255 

Table  2:  CPB06  yield  function  coefficients  corresponding  to  various  pre-strain  levels  for 
uniaxial  compression  of  clock-rolled  zirconium  along  the  through-thickness  direction. 
For  all  levels  of  effective  plastic  strain:  a  =  2  and  Li  i  =  1 .0. 


k 

L,2 

Ll3 

L22 

L23 

L33 

L44 

L55 

1-66 

0.2 

0.07 

2.883 

1.875 

1.252 

1.965 

1.204 

0.989 

1.075 

3.956 

1% 

0.23 

2.964 

1.556 

0.803 

1.476 

0.433 

1.152 

1.173 

3.313 

5% 

0.40 

2.652 

1.619 

0.813 

1.541 

0.583 

0.988 

1.127 

2.858 

25 

0.19 

3.348 

1.575 

1.039 

1.525 

0.603 

1.883 

2.081 

3.391 

35 

0.08 

3.272 

1.524 

1.120 

1.471 

0.641 

2.122 

2.312 

3.847 

Table  3:  CPB06  yield  function  coefficients  corresponding  to  rolled  tantalum;  a  =  2  and 
Ln  =  1.0. 


k 

L|2 

Ll3 

L22 

L23 

L33 

0 

-0.1911 

-0.0687 

1.0411 

0.0067 

1.1366 

39 


